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mounts. 
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The  work  documented  here  was  done  in  the  Space  and  Surface  Systems  Division  and  was 
funded  as  part  of  the  development  of  computer  programs  connected  with  the  evaluation  of  sea¬ 
mount  survey  techniques. 
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INTRODUCTION 


Except  for  certain  minor  perturbations,  the  sea  surface  is  an  equipotential  surface  in  the 
earth’s  gravity  field.  As  such  it  is  perpendicular  everywhere  to  the  gravity  vector,  which  in  the 
presence  of  a  seamount,1  2  tilts  toward  the  seamount’s  center  of  gravity.  This  tilt  in  turn  produces 
a  surface  elevation  directly  above  the  seamount. 

Being  artifacts  of  seamount  gravitation,  surface  deflections  of  this  kind  (also  called  “sea¬ 
mount  signatures”)  must  reflect  the  shape  of  and  the  mass  distribution  within  the  seamount. 
Above  the  major  seamounts,  the  maximal  geoidal  elevations  may  amount  to  several  meters  with 
corresponding  maximal  vertical  deflections  in  excess  of  10  arc  sec.  Signatures  associated  with 
major  known  seamounts  are  clearly  seen  on  the  satellite  altimetry  data  tracks.3'20  It  is  now 
widely  accepted  that  seamount  signatures  can  be  separated  from  the  noise  and  from  the  remaining 
features  on  these  data  records.  Digital  filtering  and  other  mathematical  pattern  recognition 
methods  are  the  obvious  means  to  that  end.  Yet  the  technical  implementation  of  the  principle 
must  still  be  demonstrated  by  a  successful,  published,  automatic  seamount  survey  based  on 
satellite  altimetry. 

The  present  computer  program  is  an  experimental  seamount  detector.  It  was  developed 
to  test  the  feasibility  of  exploiting  GEOS-3  and  SEASAT-A  radar  altimetry  with  the  purpose  of 
identifying  seamount  overflights,  of  locating  the  seamount  peaks  (the  seamounts  themselves  are 
several  tens  of  kilometers  wide,  typically),  and  of  “measuring”  the  submergence  depths  of  the 
peaks.  The  mission-oriented  aspect  of  this  effort  is  that  satellites  can  survey  vast  ocean  areas  much 
more  inexpensively  and  faster  than  ships  and  that  knowledge  of  the  just-mentioned  parameters 
over  large  regions  would  minimize  the  hazard  to  submarine  navigation  due  to  the  presence  of 
seamounts. 

In  essence,  this  seamount  detector  features  a  digital  filter  necessary  to  recognize  and  enhance 
the  seamount  signatures  concealed  among  the  radar  altimetry  data  (consisting  of  geoid  height 
and  vertical  deflection,  suitably  filtered),21,22  a  physical  model23  containing  the  seamount 
related  potential  theory,  and  a  mathematical  estimator  for  the  physical  seamount  parameters,  in 
particular  peak  submergence  depth.  The  latter  estimator  adjusts  the  characteristic  properties  of 
the  radar  altimetry  signature  predicted  from  the  theoretical  model  to  the  corresponding  quantities 
among  the  empirical  data.  It  is  expected  that  this  experimental  detector  algorithm  differs  from 
future,  more  advanced,  automatic  seamount  detectors  in  the  details  of  the  solution  for  the  sea¬ 
mount  model  parameters.  Specifically,  the  present  experimental  version  is  founded  on  the  as¬ 
sumption  that  the  seamount  slope  is  proportional  to  the  maximal  vertical  deflection  of  the  as¬ 
sociated  geoid  height  signature.  The  latter  relationship  had  been  devised  by  studying  some  50 
GEOS-3  radar  altimetry  data  tracks  crossing  seamounts  in  the  New  England  seamount  province. 
The  study  was  confined  to  that  part  of  the  world  because  a  dense  pattern  of  altimetry  tracks 
and  rather  reliable  bathymetry  were  both  available  there.  Once  determined  from  this  empirical 
formula,  the  seamount  slope  angle  is  retained  throughout  the  procedure.  Applying  the  above 


1 


mentioned  seamount  model,2  3  the  algorithm  then  iteratively  varies  the  remaining  seamount 
dimensions  (width  at  the  base,  peak  submergence  depth)  until  the  computed  maximal  geoid 
deviation  above  the  seamount  peak  matches  the  altimetry  signature  height  within  a  specified 
margin,  attributing  the  resulting  peak  depth  and  base  width  to  the  actual  seamount.  Because  of 
our  very  imperfect  knowledge  of  the  conditions  in  the  crust  and  upper  mantle  undenearth  sea¬ 
mounts  (problem  of  seamount  compensation),  this  procedure  is  executed  twice  per  detection, 
separately  for  the  cases  of  perfect  isostasy  and  absence  of  root.  The  actual  depth  value  is  assumed 
to  be  bracketed  by  the  depth  estimates  resulting  from  the  two  ideal  cases.  As  an  option,  the 
depth  estimation  may  additionally  be  performed  for  any  desired  degree  of  general  compensation 
with  full  freedom  in  selecting  seamount  root  height,  width,  and  density  excess  over  the  mantle. 

Computer  coding  and  computer  program  checkout  of  this  seamount  detector  have  been 
completed.  To  evaluate  its  performance,  it  is  now  being  applied  to  selected  SEASAT-A  altimetry 
tracks  in  various  oceans.  Already  it  appears  that  upon  completing  some  required  fine-tuning  of 
its  numerical  parameters,  this  computer  program  will  become  a  useful  device  for  satellite  altimetry 
data  analysis. 

One  more  detail  should  finally  be  mentioned.  It  was  originally  intended  to  write  a  seamount 
detector  capable  of  fully  automatic  operation.  But  the  numerical  program  checkout  as  well  as 
the  current  evaluation  using  actual  altimetry  indicate  that  any  operational  success  for  seamount 
diagnostic  would  be  considerably  enhanced  by  human  intervention.  Parallel  to  the  computer 
coding  effort  for  the  automatic  detector,  a  computer  graphics  system  was  thus  developed,  which 
permits  the  visual  scanning  of  altimetry  data  over  large  ocean  areas.  It  generates  charts  of  suitably 
filtered  and  edited  altimetry  that  may  be  overlaid  on  bathymetric  charts  where  such  are  available. 
Use  of  these  chart  data  appears  to  greatly  increase  the  effectiveness  of  the  experimental  automatic 
seamount  detector.  In  fact,  they  may  turn  out  to  be  indispensable  for  the  success  of  any  seamount 
survey  now  at  least  and  in  the  near  future. 


EXPERIMENTAL  COMPUTER  ALGORITHM  FOR 
SEAMOUNT  MODEL  PARAMETER  ESTIMATION 

SATELLITE  RADAR  ALTIMETRY  DATA  FORMAT 

The  satellite  altimetry  data  to  be  processed  by  the  seamount  detector  are  available  from  per¬ 
manent  storage  in  the  formats  specified  in  Tables  1  through  4  (G.  B.  West,  unpublished  data, 
1978).  Select  the  type  of  altimetry  desired,  GEOS-3  or  SEASAT-A,  and  extract  from  the  header 
of  the  corresponding  data  file  the  time  interval.  At,  in  seconds  and  the  subsatellite  velocity,  vs, 
in  km/sec.  From  each  segment  of  the  altimetry  data  record,  obtain  for  each  data  point  (time 
instant,  tif  in  seconds)  the  filered*  geoid  height,  Ni%  =  N(tj),  in  meters,  the  deflection  of  the 
vertical,  8X  =  5(tj),  in  arc  seconds,  the  latitude,  ^  =  <p(t|),  in  radians  and  the  longitude  East, 
Xj  =  X(tj),  in  radians.  Convert  ^  and  \  to  decimal  fractions  of  degrees. 


♦The  reader  is  urged  to  keep  in  mind  that,  as  said  in  the  INTRODUCTION,  the  altimetry  data  serving  as  input  to  this  seamount 
detector  are  not  raw  altimetry  data  tracks  but  result  themselves  from  a  filtering  process  explained  in  References  21  and  22.  Here 
and  in  the  following  we  shall  refer  to  Nj  and  6j  as  “filtered”  when  wishing  to  emphasize  their  origin,  by  Kalman  smoothing,  from 
the  raw  satellite  radar  altimetry.  With  reference  to  the  present  seamount  detector,  Nj  and  6$  will  usually  be  considered  “unfUtered” 
and  named  so  as  they  are  input  data  to  the  matched  filter. 
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Table  1.  G EOS-3  Satellite  Radar  Altimetry  Data  Filtered  Along  Track  Geoid  Heights 
and  Vertical  Deflections  (Sinary  File  Format,  Header  Record) 


Word  Number 

in  Record 

Format 

Description 

1 

I 

Filter  word  2  «  Filtered,  3  =  linear  fit 

2 

I 

Kalman-Wiener  control  -  0  -  Wiener  Filter  output 

1  -  Kalman  Filter  output 

3 

I 

Number  of  data  points  in  track 

4 

l 

Measurement  type  (40  =  G  or  41  =  H) 

5 

1 

Time  system  indicator 

Beginning  time  of  data 

6 

F 

Year 

7 

F 

Day  of  year 

8 

F 

Time  of  day  (seconds  from  beginning  of  day) 

9 

F 

Time  interval 

10 

F 

Subsatellite  velocity  (km/sec) 

Ending  time  of  data 

11 

F 

Year 

12 

F 

Day  of  year 

13 

F 

Time  of  day  (sec  from  beginning  of  day) 

14 

A 

Date  of  raw  data  tape 

15 

A 

Date  of  reduction 

16 

F 

Geoid  height  mean  value  over  data  track 

17 

F 

Standard  deviation  with  respect  to  the  mean  for  geoid  heights  (M) 

18 

F 

Deflection  of  the  vertical  mean  value  over  data  track 

19 

F 

Standard  deviation  with  respect  to  the  mean  for  deflections  of  the  vertical 

20 

I 

Processing  report 

1-9  Ex,  ellipsoid.  Filtering  parameter,  etc. 

21 

I 

Mode 

22 

I 

Data  status  word 

23 

1 

Rev  number 

24 

I 

Total  number  of  points  on  File  (entire  pass) 

25 

I 

Number  of  segments  in  pass 

26 

1 

Station  number 

27 

1 

Tide  indicator  0  •  Before  tide  correction 

I  -  After  tide  correction 

28 

F 

Autocorrelation  distance 

29 

F 

Standard  deviation  of  the  noise 

30 

F 

Standard  deviation  of  the  vertical  deflection 

31 

F 

Geoid  height  bias  estimate  (m) 

32 

F 

Vertical  deflection  bias  estimate  (arc  sec) 

33 

F 

Vertical  deflection  uncertainty  estimate  (maximum) 

34 

F 

Vertical  deflection  uncertainty  estimate  (steady-state) 
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Table  2.  GEOS-3  Satellite  Radar  Altimetry  Data  Filtered  Along  Track  Geoid  Heights 
and  Vertical  Deflections  (Binary  File  Format,  Packed  Data  Record) 


Word  Number 

in  Record 

Format 

I 

I 

i 

1 

Description 

Filtered  geoid  height  (in)  (FGH)+ 
raw  geoid  height  (in)  ( RGH )  + 
tide  correction  (in)  (TO  +  geoid  height 
uncertainty  estimate  (ni)  (GHUFj 

Packed: 

KGH+500)  x  100]  x  10"  +  KRGH+500)  x  10]  x  10s  + 
[(TC+5)  x  100]  x  102  +  IGHUE  x  100] 

Detleetion  of  vertical  in  arc  sec  (DV)  + 
latitude  in  rad  (xLAT )  + 
longitude  in  rad  (xLON) 

Packed: 

KDV+500)  x  10]  x  10'°  +  [(xLAT+5)  x  104]  x  10s 
lx  LON  x  104] 


Table  3.  SEASAT-A  Filtered  Geoid  Height  Data  File  (Type  FGD) 
Header  Record  (Common  Header) 


Approximate  Range 
Word  Type  of  Significance 


Description 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


I 

l 

I 

l 

R 

R 

A 

1 

R 

R 

R 

R 

R 

R 

R 

R 

K 

R 

R 


xxxx 

xxxxx 

XX 

XXX 

XXXXX.  XXX 

XXXXXX 

YYDDD 
XXX 
XXXX 
XX. XX 
XX  XX 
XX  XX 
X  XX 
-XXX 

X  XXX 

x.xx 

XX. X 
t.XXX 


Number  of  points  in  file 
Rev  number 

Starting  year  of  segment 
Starting  day  of  segment 
Starting  seconds  of  segment  (sec) 

Time  interval  tsec) 

Year  and  Julian  day  (right  adjusted) 

Altimeter  mode 
Autocorrelation  distance  (km) 

Standard  deviation  of  data  (m) 

Standard  deviation  of  geoid  heights  (m) 

Standard  deviation  of  vertical  deflections  (arc  see) 
RMS  of  fihered-raw  differences  (in) 

Maximum  vertical  detleetion  (arc sec) 

Average  velocity  (km/sec) 

Antenna  distance  from  satellite  center  of  gravity  (mi 
Radar  instrument  delay  distance  equivalent  (cm) 
Time  correction  <see> 

Spare 
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Table  4.  SEASAT-A  Filtered  Geoid  Height  Data  File  (Type  FGD)  Data  Record  (Common  Data) 


Approximate  Range 


Word 

Type 

of  Significance 

1 

R 

+X.XXXXX 

■> 

R 

x.xxxxx 

3 

R 

+XXX.XX 

4 

R 

±xxx.xx 

5 

R 

±xxx.xx 

6 

R 

x.xx 

7 

R 

xx.xx 

8 

R 

x.xx 

9 

R 

XXX 

10 

1 

I,  --1,4 

11 

R 

xx.x 

12 

R 

xx.x 

13 

R 

XX 

14 

R 

XX 

15 

R 

±.xx 

16 

R 

x.xxx 

17 

R 

xx.x 

18 

R 

XXX 

19 

R 

xxx.x 

20 

R 

xx.x 

21 

R 

±xx.xx 

7*> 

R 

x.xx 

23 

R 

XX 

24 

R 

XX 

25 

R 

±x.xx 

26 

R 

xx.x 

27 

R 

xx.x 

28 

R 

xxx.x 

29 

R 

30 

R 

Description 

Latitude  (rad) 

Longitude  (rad) 

Geoid  height  (m) 

Vertical  deflection  (arc  sec) 

Raw  geoid  height  (m) 

Geoid  height  confidence  bound  (m) 

Vertical  deflection  confidence  bound  (arc  sec) 
Orbit  uncertainty  in  geoid  heights  (m) 

Orbit  uncertainty  in  vertical  deflections  (arc  sec) 
Flag  word 

Significant  wave  height  (m) 
o  -  SWH  (m) 

Automatic  gain  control  (dB) 
a  -  AGC  (dB) 

Tilt/SWH  correction  (m) 

Tilt  (rad) 

Ionospheric  correction  (cm) 

Atmospheric  pressure  (mb) 

Dry  tropospheric  correction  (cm) 

Wet  tropospheric  correction  (cm) 

Tide  (m) 

Barotropic  correction  (m) 

Wind  speed  (kn) 

Wind  direction  (rad) 

Sea  state  (SWH)  correction  (m) 

Alternate  wet  tropospheric  corection  (cm) 

Rain  rate  (mm/hr) 

Steric  correction  (cm) 

Spare 

Spare 
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MATCHED  FILTER 


The  matched  filter  (B.  Zondek,  unpublished  data,  1980)  is  an  optimized  high-pass  filter. 
To  apply  it  to  SEASAT-A  data,  the  control  card  deck  specified  in  Table  5  is  requircu.  Also,  six 
data  cards,  “A”  through  “F,”  must  follow  the  control  car  1  deck.  The  latter  contains  two  SCOPE 
cards  selecting  the  Revolution  Number.  By  way  of  an  example  indicated  by  pencil  arrows,  a 
specific  Device  Set  Number  and  a  File  Name  are  inserted.  The  six  data  cards  determine  the  re¬ 
quired  calculations.  They  may  be  manipulated  as  follows. 

Card  A  (selecting  the  filter  type): 

Col.  6-10:  Enter  “3”  into  Col.  10. 

Replace  this  by  a  "0,”  if  printout  of  unfiltered  data  plus  counting  numbers  wanted. 
Col.  11  -  13:  Starting  time  Tj  in  sec. 

Col.  31-50:  Terminating  time  T2  in  sec. 

Card  B  (Use  or  omission  of  starting  and  terminating  time): 

Col.  5:  Enter  “0,”  if  Tt  and  T2  to  be  ignored. 

Enteral,”  if Ti  andT2  to  be  used. 

C  ard  C  (Choice  of  printout): 

Col.  5:  Enter  “0”  for  extended  printout  (14  digits  for  results  and  intermediate  results). 

This  is  a  program  checkout  mode.  Enter  “1”  to  facilitate  printing,  plus  writing  on 
a  tape,  of  restricted  output.  This  is  the  standard  output  format. 

Enter  “2,”  if  desiring  printing  (only)  of  restricted  output  with  e  choice. 

Enter  “3”  to  print  and  write  on  tape  restricted  output  with  e  choice. 

Card  D  (Specifying  number  of  filter  tuning  constants): 

Col.  1-5:  Enter  number  of  filter  tuning  constants  p.  For  each  p,  the  computer  program 

will  process  and  print  out  the  original  SEASAT-A  data  track. 

Card  E  (Specifying  the  filter  tuning  constants): 

Col.  1  -  10: 

Col.  11-20:  Card  may  contain  up  to 
Col.  21-30:1  eight  different  p  values  in 
etc.  F  10.9  format. 

Col.  71  -  80:  j 

Card  F  (e  card): 

Given  a  value  for  e,  this  card  will  prevent  printing  of  a  particular  time  line  (geoid  height  and 
deflection  of  vertical  versus  time  in  sec),  if  151  <  e.  Reasonable  values  for  e  are  0.1  <  e  <  10. 
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Adopt  the  following  terminology:  N.  are  the  “unfiltered  geoid  height”  values,  and  are 
the  “unfiltered  vertical  deflection”  values.  Application  of  the  matched  filter  to  N  and  6  will 

A  A  A  I  a  1 

result  in  the  NT  and  5.  data  tracks.  N.  are  the  “filtered  geoid  heights”  (in  m)  and  6  are  the 
“filtered  vertical  deflections”  (in  arc  sec).  The  index  “i”  is  related  to  the  time  value  t  associ¬ 
ated  with  the  data  values. 

Note  that,  because  of  the  nature  of  the  filter  algorithm,  no  filtered  data  values  vill  be  as¬ 
sociated  with  the  first  few  and  the  last  few  unfiltered  data  points.  Thus,  the  first  “da*.,  window” 
(see  below  for  definition)  on  a  data  track  shall  start  with  the  first  available  filtered  geoid  height 
and  deflection  of  the  vertical  and  the  last  data  window  shall  end  with  the  last  available  pair  of 
filtered  data  values. 

Application  of  the  matched  filter  to  GEOS-3  satellite  altimetry  is  omitted  from  this  report. 

Table  5.  Control  Card  Deck  and  Data  Card  Listing  for  Matched  Filter 


USER  NUMBER  (T400,  P3,  PL-1( 

CHARGE  CODE,  PASS  WORD. 

ATTACH,  OLDB1N,  BHPSABIN02,  ID=N60,MR=1 . 
BEGIN,  MOUNT,  VSN=|NUA022. 1 - 


ATTACH,  TAPE  1,  FGD01 166,  SN=NUA022,1D=N14,MR=1. 
ATTACH,SSL1B,1D=NV5. 


ATTACH.  SYSLIB. 
LIBRARY,  SYSL1B,  SSLIB. 
LDSET  (PRESET=INDEF) 
OLDBIN. 

EXIT. 

DMP  (0,2(02,421) 

7/8/9 

DATACARD  "A” 
DATACARD  “B” 
DATACARD  “C” 
DATACARD  “D” 
DATACARD  “E” 
DATACARD  “F” 

6/7/S/9 


File  Name 


Device  Set  Number 


ROUGHNESS  DETECTOR 


Input  Data 

A 


Ni 

Data  track  (filtered  geoid  height) 

m 

A 

6. 

1 

Data  track  (filtered  vertical  deflection) 

arc  sec 

e 

Parameter  of  matched  filter 

arc  sec 

At 

Data  step  width.  At  =  t.  -  t._  x 

sec 

vs 

Subsatellite  velocity 

km/sec 

L 

Data  window  width  (specified  by  user) 

km 

deln 

Threshold  parameter  (specified  by  user) 

m 

Algorithm 

A  A 

A 

1 .  Consider  the 

)  filtered  data  tracks  N.  and  6  .  Assign  the  value  zero 

to  all  N.  that  are 

negative.  Specify  data  track  44 windows.”  The  number  of  the  data  points  within  each 

window  is 

n  =  INT  f  \ 

+ 1 

(301) 

\vs  At/ 

Start  the  first  window  at  the  first  data  point  not  omitted  by  the  filter, 

tfc,  and  calculate 

/i  «  / 

A 

v 

sn‘=V"S( 

V 

,+i) 

(302) 

Do  the  same  for  the  next  window,  NR+1 ,  NR  +  2, ...  ,  NR+n,  obtaining 

h  Vi 

/  * 

v 

(N‘ 

-,+i) 

(303) 
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and  for  each  following  window,  obtaining 


SN 


k  +  v 


(304) 


until  the  data  track  N.  is  exhausted. 

2.  From  now  on,  consider  only  those  windows  for  which  SNV  >  deln.  From  the  remaining 
windows,  select  those  for  which  SNV  is  larger  than  the  SN  values  for  the  neighboring  windows. 
(In  case  the  set  SNV  =  fct(v)  does  not  have  a  single  peak  but  a  “flat  top,”  select  that  window  that 
is  associated  with  the  middle  of  the  flat  top).  From  here  on,  the  term  “window”  shall  refer  to  a 
member  of  the  just  defined  set  of  maximum-SNv-windows. 


Output  Data 


For  each  of  these  windows  (characterized  by  maximum  SNV  value),  store  and  print  out 

+ 1  +  n  ~  1 

A  A  A 

M  M+ 1  M+n“l 

A  A  A 

+  1  .  ^M  +  n-1 

SNV 

Nm  Nm+1  .  NM+n_j 

^M+l  .  ^>l  +  n-l 


Additional  Details 

1 .  To  avoid  disrupting  the  continuity  of  the  data  record  across  any  seamount  signature, 
assign  to  the  filter  parameter  epsilon  the  value  zero. 

2.  This  means  that  the  window  search  will  now  be  conducted,  without  interruption,  from 
the  beginning  to  the  end  of  any  continuous  portion  of  a  particular  data  track. 

3.  Eliminate  all  data  track  portions  that  after  application  of  the  matched  filter  have  a 
number  points  which  is  less  than  1 .5  to  2.0  times  the  number  n  of  data  points  within 
the  window,  Equation  (301). 
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SEAMOUNT  LOCATOR 


Subroutine  “EXTR” 

Known  are  three  data  points,  i  =  1,  2,  3,  on  a  data  track  y;  =  y(Xj).  Assume  that  these  points 
bracket  a  maximum  or  minimum  of  the  yf  values.  The  abscissa  value,  x  ,  associated  with  the 
extremum  of  y.  is  then 


(401) 


y2  -  yi  (x2-x,)2 

bd  =  (402) 

y3  -  yi  (x3-x,)2 


x2  -  x,  y2  -  y1 

cd  -  (403) 

x3-x,  y3  -  yi 


Subroutine  “INTERPOL-1” 

Known  are  three  data  points,  i  =  1,  2,  3,  on  a  data  track  y.  -  y(Xj).  Find,  for  a  specified 
abscissa  value,  xp,  the  ordinate  value  yp  =  y(xp). 


y  =  a  +  b  (x„  - 

Jp  P 

Xj)  +  c(xp  -  x,)2 

(404) 

a  =  y, 

(405) 

b.i 

y2  -  yi 

(x2  -  x,)2 

(406) 

d 

ya  -  yi 

(X3  -  X,)2 

1 

X2  -  x, 

y2  -  yi 

(407) 

d 

x3  - 

y3  -  yi 

d  =  (x2  -  x, )  (x3  -  x, )  (x3  -  x2)  (408) 
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Subroutine  'iNTERPOL-2” 


Known  are  tour  data  points,  i  =  l,  2,  3,  4,  on  a  data  track  yj  =  y(x/).  Calculate,  lor  a 
specified  abscissa  value,  xp ,  the  ordinate  value  yp  =  y(xp  ). 


yp  =  c  +  f(*P 

-  x,)  +  g(xp 

-  Xj  )2  +  h(xp 

-x,)3  (409) 

e 

(410) 

i 

y 2  -  yi 

(X 2  -  X, 

)2 

(x2  -  X, 

)3 

ii 

t>  l- 

yj  -  yi 

(X3  -  X, 

)2 

(x3  -  X, 

)3 

(41 1 ) 

y4  -  yi 

(x4  -  X, 

(x4  -  X, 

)3 

1 

x2  -  X, 

y3  -  yi 

(x2  -  X, 

)3 

g" 

x3  -  X, 

y3  -  yi 

(x3  -  X, 

)3 

(412) 

x4  -  X, 

y4  -  yi 

(x4  -  Xj 

)3 

1  x2  -  Xj 

(x2  -  X, 

)2 

y3  -  yi 

1 

1 

h  =  - 
A 

X3  "  X1 

(x3  -  X, 

)2 

y3  -  yi 

(413) 

x4  ”  X1 

(x4  -  X, 

)2 

y4  -  yx 

i 

X2  -  X! 

(x2  -  X, 

)2 

(x2  -  X, 

)3 

4= 

x3  “  xl 

(X3  -  X, 

)2 

(x3  -  X, 

)3 

(414) 

x4  -  X, 

(x4  -  X, 

)2 

(x4  -  X, 

)3 

Note  Concerning  the  Subroutines 

Optionally,  the  computer  programmer  may  implement  the  above  subroutines  by  any  of  the 
classical  interpolation  schemes.  If  desired,  the  polynomial  coefficients  can  be  evaluated  by  a  least- 
squares  fit  method.  Irrespective  of  the  method  chosen,  it  must  however  be  kept  in  mind  that 
rather  frequently  along  the  altimetry  track,  the  argument  values  (time  or  distance  -  depending  on 
the  interpretation  of  x.)  will  be  large  compared  to  the  time  required  for  a  seamount  overflight  or 
to  the  seamount  width,  implying  x.  >  (xk  -  x.).  To  avoid  loss  of  significant  digits  during  the 
evaluation  of  the  determinants,  it  is  advisable  to  first  transform  the  x{  values  to  corresponding 
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values  so  that  the  origin  of  the  x  system  coincides  with  ,  x.  =  (xi  -  Xj ),  implying  Xj  =  0. 
The  particular  interpolation  method  chosen  may  then  be  executed.  Subsequently,  the  results 
must  be  transformed  back  into  the  x/y  coordinate  system.  Note  that  this  precaution  is  unnecessary, 
if  the  subroutines  are  applied  in  the  form  specified  above. 


input  for  Seamount  Locator  Algorithm 


For  each  data  window,  have  available 


w 

Unfiltered  geoid  height 

m 

w 

Unfiltered  vertical  deflection 

arc  sec 

t 

1 

Time  argument 

sec 

xi 

Longitude  East  from  altimetry  file 

deg 

Latitude  from  altimetry  file 

deg 

Seamount  Locator  Algorithm 

For  each  of  the  specified  data  windows,  now  consider  the  unfiltered  geoid  heights,  Nt,  and  the 
unfiltered  deflections  of  the  vertical,  6^  Also  keep  in  mind  that  it  is  possible  to  associate  with 
the  discrete  number  sets  N.(tj)  and  the  continuous  functions  N(t)  and  <5(t)  representing  the 

geoid  height  and  vertical  deflection  in  the  “real  world”  or  simply  being  assumed  to  be  suitable 
functions  adapted  to  N.(t.)  and  by  least-squares  fit  or  other  meaningful  methods.  For  the 
present  purpose,  it  is  sufficient  to  postulate  that  N(t)  and  5(t)  may  be  realized  mathematically  if 
needed.  Perform  the  following  calculations. 

1.  When  traversing  the  window  from  the  left  to  the  right  (proceeding  from  lower  values  of 
t.  to  higher  ones),  expect  to  encounter  negative  5i  values  and,  among  them,  an  absolute 
minimum  value.  Consider  the  data  point  5.(t.)  associated  with  this  minimum  value  plus 
the  data  point  to  the  left  and  the  data  point  to  the  right  of  it.  Apply  to  these  three  data 
points  the  subroutine  “EXTR”  to  find  the  abscissa  value  tA  belonging  to  the  minimum  of 
the  curve  5(t)  associated  with  the  just  specified  three  data  points. 


In  the  same  manner,  find  the  time  value  tB  associated  with  the  absolute  maximum  value 
of  6:(t.)  that  will  occur  further  down  the  data  track  but  still  within  the  data  window. 
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3.  Within  the  window,  the  data  track  N4(t)  may  be  expected  to  have  an  absolute  maximum. 
Find  the  data  value  representing  this  maximum.  Also  find  the  neighboring  data  points 
to  the  left  and  to  the  right  of  this  maximum.  Apply  subroutine  “EXTR”  to  find  the  time 
value  t$M  of  the  maximum  of  N(t)  related  to  the  just-defined  three  data  points. 

4.  Now  apply  subroutine  “INTERPOL-1”  to  tSM  and  the  three  data  points  on  the  N. 
track  from  which  tSM  was  determined  in  step  3,  obtaining  N(tSM )  =  NSM  . 

5.  Next,  determine  the  four  data  points  on  the  Ns  track  that  are  nearest  to  the  time  value 
tA .  Apply  to  them  subroutine  “INTERPOL-2”  to  evaluate  the  N  value  associated  with 
tA,N(tA)  =  NA. 

6.  Further,  proceed  as  in  step  5  to  find  N(tB )  =  NB . 

7.  Proceed  as  in  step  5  to  find  6(tA)  =  6a. 

8.  Proceed  as  in  step  6  to  find  8(tfi )  *  5B . 

9.  Consider  now  the  data  tracks  and  6.  in  the  neighborhood  of  tSM ,  in  particular,  the 
time  value  t.  =  tjM  just  to  the  left  of  (prior  to)  tSM  and  the  time  value  =  t*M  to  the 
right  of  (following)  tSM.  Note  that  t”M  <  tgM  and  t*M  >  tSM .  In  the  radar  files,  t~M 
and  t J  M  are  associated  with  the  values  of  longitude  East  and  latitude  on  the  subsatellite 
track, 

*SM  ''  XSM  3nd  *SM 

lSM  ~  XSM  and*SM 

Evaluate  the  longitude  East  and  the  latitude  of  the  “signature  peak,” 

XSM  =XSM  (415) 

t+  -  t" 
lSM  lSM 

M  “  ^SM 

^SM=^SM+_ ~  ~  (tSM_tSM^  (416) 

(SM  "  fSM 
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Wk^'l 


10.  Note  that  above  the  data  were  postulated  to  conform  to  the  typical  seamount  altimetry 
pattern  illustrated  in  Figure  1 .  If  the  data  fail,  altogether  or  in  any  essential  detail  to 
match  this  pattern,  reject  the  data  window. 


Seamount  Locator  Output 


For  each  data  window  successfully  processed,  store  and  list  the  following. 


^B  ^B 

t  N 

lSM  nSM 


M  M 


Seamount  longitude  and.  latitude  are  to  be  specified  in  decimal  degrees  as  well  as  in  degree  - 
minute  -  second  format. 


CALCULATION  OF  INITIAL  SEAMOUNT  PARAMETERS 

From  the  data  output  of  the  preceding  computer  program  segments,  calculate  the  following 
quantities  for  each  data  window  successfully  processed. 


W„-2B„  »*•«,-«*> 
W 

B.~2 


Nt“l-7[NSM4<N»+N*)] 


a 


7T 


A  180  3600 


Ofn  + 


7T  5B 


B  180  3600 


(501) 

(502) 

(503) 

(504) 

(505) 
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VERTICAL  DEFLECTION  (arc  sec)  GEOID  HEIGHT  (m) 


**  (1.5E  +  05)  aA 


(506) 


*>B  *  (1.5E  +  05) 


(507) 


Wy  is  the  initial  estimate  of  “seamount  width  at  the  base”  in  km.  BQ  is  the  “seamount  half  width 
at  the  base”  in  km.  Nc  is  the  geoid  elevation  “observed”  at  the  seamount  center,  in  m,  to  be 
matched  by  the  calculated  geoid  elevation  during  parameter  adjustment.  y?A  and  are  the 
estimated  seamount  slope  angles  in  deg.  <*A  and  aB  are  the  maximal  signature  slope  angles  in  rad. 


OCEAN  DEPTH  ALGORITHM 

For  each  data  window  successfully  processed,  the  ocean  depth,  D,  at  the  seamount  perimeter 
is  required.  Suitable  ocean  depth  data24  are  available  as  averages  over  one-by-one  degree  surface 
area  elements,  as  a  permanent  data  file  on  the  Dahlgren  computer  system.  If  necessary,  convert 
D  to  meters. 


ESTIMATION  OF  PEAK  DEPTH 


Subroutine  “DN”  -  Input  Data 


D 

Ocean  depth  at  seamount  perimeter 

m 

ds  or  d 

Depth  of  seamount  peak 

m 

Seamount  slope  angle 

deg 

Hs 

Seamount  height 

m 

T 

Crustal  thickness 

nr 

Br 

Half  width  of  seamount  root  as  base 

m 

hr 

Depth,  below  crust,  of  seamount  root 

m 

Ps 

Seamount  density 

gr/m3 

Pw 

Water  density 

gr/m3 

Pr 

Root  density 

gr/m3 

Pm 

Mantle  density 

gr/m3 

G/g 

Ratio  of  Newton's  constant  to  surface  gravity 

m2/gr 

If  not  otherwise  specified,  assume 


16 


Ps  =  2.60  E+06  gr/m3 

pw  =  1.03  E+06  gr/m3 

pR  =  2.95  E+06  gr/m3 

pM  =  3.40  E+06  gr/m3 

G/g  =  0.68024  E-14  m2/gr 
Subroutine  “DN”  -  Algorithm 

Reference  is  made  to  Figure  2.  Perform  the  following  calculations. 

as  =  tantps  (701) 

jSs  =  d/Hs  (702) 

aR  =  Hr/br  (703) 

0r  =  (D  +  T+Hr)/Hr  (704) 

evaluate  FU(as, /3S)  (Subroutine  “FU”)  (705) 

evaluate  FI  (aR ,  0R)  (Subroutine  “FI”)  (706) 

DNS  =  (p$  -  pw  )  H|  j  FU(as ,  ps )  (707) 

DNR  =  (pM  -  pR)HR  j  FI(aR,0R)  (708) 

DN  =  DNS  -  DNR  (709) 

Subroutine  “DN”  -  Output 

DNS  in  m 
DNR  in  m 
DN  in  m 


17 


Subroutine  “FU” 


FU(a,  p)  - 


=  n t 


-a2  )32 
2(1  +  a2 ) 


1  +  a2(l  +  p) 
2a(  1  +  a2 ) 


\/l  +  a2  (1  +  P)2 


a  p2  /  a  P  ( yj  1  +  a2  -  a) 

2(l+a2)3/2  y^/T+a2”  ^/1+ce2  ( 1+/3)2  -  1  -  a2  (1  +0) 

1 


Subroutine  “FI” 

FI(«,  P)  = 


a2  ff2  1  -  a2(/3  -  1) 
2(1  +  a2 )  2a(l  +  a2 ) 


Vl  +  a2(0-  l)2 


+  ^  / _ ag(  %/l+a2  +  a) _ 

2(l+“2)3/2  \\/l+a2  \/  1+or2 (/3-  1  )2  -  l-a203-  1) 
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AXIS  OF  ROTATIONAL 
SYMMETRY 
I 


SEA  SURFACE 


Figure  2.  Geometry  of  Seamount  Model 
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Estimation  of  Peak  Depth  for  Isostaticallv  Compensated  Seamount-Input  Data 


(i/g\  Ratio  of  Newton’s  constant  to  surface  gravity  (user  specified)  nr  gr 


Ps  I  Seamount  density  (user  specified)  gr/m3 

pvv  \  Water  density  (user  specified)  gr/m3 

pR  V  Root  density  ( user  specified)  gr/m3 

pM  /  Mantle  density  (user  specified)  gr/m3 

D  Ocean  depth  at  seamount  perimeter  m 

T  C  rustal  thickness  (user  specified)  m 

,  Seamount  slope  angles  deg 

V 

Initial  estimate  of  seamount  width  at  base  km 

Nc  “Observed"  geoid  elevation  above  seamount  peak  m 


If  not  otherwise  specified,  assume  for  G/g  and  the  four  densities,  the  nominal  values  listed  under 
“Subroutine  “DN"  -  -  Input  Data." 


Estimation  of  Peak  Depth  for  Isostatically  Compensated 
Seamount  -  -  Removal  of  Ill-Conditioned  Cases 


Reference  is  made  to  Figure  2.  Execute  the  following  test.  This  test  assumes  that  the 
seamount  is  as  large  as  the  estimated  slope  angle  *ps  and  the  local  ocean  depth,  D  permit  (peak  is 

near  surface . d*  -  10  ml,  thus  producing  the  largest  maximal  geoid  elevation  DN*  possible 

for  the  particular  and  D  values.  If  DN*  is  smaller  than  the  observed  maximal  geoid  elevation 
Nc ,  the  next  program  segment  will  be  unable  to  perform  the  depth  estimation  unless  certain 
corrective  measures  are  taken  as  prescribed  below. 
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^A  + 


tig  =  1 0  in  (713) 

*  D-  10  m 

LL  = -  (714) 

b  tain  i£s 

H*  =  D  -  10  m  (715) 

b‘r  =  b; 

H*  =  Hs  (Ps^fwj  (717) 

DN*  =  DN(D,d*,Y^H*,T.  dr,Hr)  (718) 


DN*  >  Nc  ?  (7i9) 

If  yes,  proceed  witn  the  next  computer  program  segment.  Otherwise,  print  out  “CAUTION: 
ILL-CONDITIONED  CASE,"  perform  the  assignment  Nc  =  DN*.  and  then  proceed  with  the 
next  computer  program  segment. 


Estimation  of  Peak  Depth  for  Isostatically  Compensated 
Seamount  -  -  Algorithm 


Reference  is  made  to  Figure  2.  Perform  the  following  claculations. 


^A  +*B 

=  2 

o 

r  | 
r 

Bso  =  500 W„ 

(721) 

do  =  D“  dso  tan*s 

<■’::) 

d0  <  0  '• 

(723) 

If  yes,  execute  Equations  (724)  and  (725)  and  continue.  Otherwise,  ignore  these  two 
equations  and  continue. 
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MAXIMAL  GEOID  HEIGHT 


i 


Figure  3.  Linear  Scheme  for  Estimation  of  Seamount  Half  Width  (at  Basel 
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(740) 


r 


"r2 


'^S2  ^ 


DN,  =  DN(D,d2,^s,HS2,T,iiR2,iIR2) 

and  use  B,, , ,  DN, ,  B,,, ,  and  DN,  to  start  the  iterative  routine 

B 


BS  i  +  2 


=  B  .  +  -LiLl — —  (N„  -  DN  ) 
s  1  DNj+1  -  DN;  c 


^i+2  ^  -  Bs  j+2  tan  < ps 

dj+  2  <  0? 


(741) 

(742) 

(743) 

(744) 


If  yes,  execute  Equations  (745)  and  (746)  and  continue.  Otherwise,  ignore  these  two 
equations  and  continue. 


d,+  2  = 

10  m 

(745) 

BS  i+2 

D-  10  m 

tan 

(746) 

HS  i  +  2 

=  Bs  i+2  tany>s 

(747) 

BR  i+2 

BS  i+2 

(748) 

HR  i+2 

i+2  ^ 

(749) 

DN,+  2 

-  DN(D.  ^i+2’^s’  i+2’  BR  i+2’  ^R  i+2^ 

(750) 

Ci+2  = 

lNc  -  dn,.2| 

(751) 

ei+2  < 

0.0001  m  ? 

(752) 

If  yes.  stop  and  call  the  associated  dj+2  “d..“  Exit  from  the  iterative  routine.  Otherwise, 
reenter  the  iterative  routine. 


End  of  the  iterative  routine. 


The  estimated  depth  of  the  seamount  peak  is 
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disoSt.  comp.  =  BSj  tan*S  =  dj 


where 


®Sj  **Si+2 


Estima tion  of  Peak  Depth  for  the  Generally  Compensated 
Seamount  -  -  input  Data 


(753 1 


( 754) 


Ratio  to  Newton’s  constant  to  surface  gravity  (user  specified) 
Seamount  density  (user  specified) 

Water  density  (user  specified) 

Root  density  (user  specified) 

Mantle  density  (user  specified) 

D  Ocean  depth  at  seamount  perimeter 

T  Crustal  thickness  (user  specified) 

*a) 

\  Seamount  slope  angles 


m2/gr 

gr/m3 

gr/m3 

gr/m3 

gr/m3 

m 

m 

deg 


WQ  Initial  estimate  of  seamount  width  at  base  km 

Nc  “Observed”  geoid  elevation  above  seamount  peak  m 

sk  Geometry  factor  (user  specified) 

Hr  “Height”  of  seamount  root  (user  specified)  m 

If  not  otherwise  specified,  assume  for  G/g  and  the  four  densities  the  nominal  values  listed  under 
“Subroutine  DN  -  -  Input  Data.” 
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Estimation  of  Peak  Depth  for  Generally  Compensated  Seamount  -  -  Removal 
of  Ill-Conditioned  Cases 

Reference  is  made  to  Figure  2.  Execute  the  following  test.  This  test  assumes  that  the 
seamount  is  as  large  as  is  compatible  with  the  estimated  seamount  slope  angle  <ps  and  the  local 

ocean  depth  D  (peak  is  near  surface . d*  =  10  m),  thus  producing  the  largest  maximal  geoid 

elevation  DN*  possible  for  the  particular  <£s  and  D  values.  If  DN*  is  smaller  than  the  observed 
maximal  geoid  elevation  Nc,  the  next  program  segment  will  be  unable  to  perform  the  depth 
estimation  unless  certain  corrective  measures  are  taken  as  prescribed  below. 


2 

755) 

d*  =  10  m 

(756) 

*  D  -  10  m 

b;  = 

(757) 

"  tan 

H*  =  D  -  10  m 

(758) 

K  =  b; 

(759) 

Hr  “  Hr 

(760) 

DN*  =  DN(D,  B^.H*) 

(761) 

DN*  >  Nc  ? 

(762) 

If  yes,  proceed  with  the  next  computer  program  segment.  Otherwise,  print  out:  “CAUTION: 
ILL-CONDITIONED  CASE,”  perform  the  assignment  Nc  =  DN*,  and  then  proceed  with  the 
next  computer  program  segment 


Estimation  of  Peak  Depth  for  Generally  Compensated  Seamounts  -  -  Algorithm 


Reference  is  made  to  Figure  2.  Perform  the  following  calculations. 
*A  +  *B 

*s  = — ; — 


(763) 
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(764) 


“so  =  500  Wo 

d0  =  D  -  BS(J  tan*s 

do  <  0  •> 

If  yes,  execute  Equations  (767)  and  (768)  and  continue.  Otherwise,  ignore 
equations  and  continue. 

dQ  =  1 0  m 

_  D-  10m 
so  tan  i^s 

^so  “so  tan  ^s 
“ro  “so 

DN0  =  DN(D.d0,^s,HS0,T,BR0,HR) 

Also,  calculate 

“si  ~  “so 
d,  =  D  -  Bs,  tan <ps 

HS 1  =  “si  tan*s 
“r  I  “s  I 

DN,  =  DN(D,  d, ,  ^>s ,  H$ , ,  T,  BR  , .  HR ) 
and  estimate  the  seamount  half  width  according  to  Figure  3. 


“s l  “so 

B.,  =  BQn  + - -  <N..  -  DN,. ) 

S2  sn  DN  -  DN  (  0 


(765) 

(766) 
these  two 

(767) 

(768) 

(769) 

(770) 

(771) 

(772) 

(773) 

(774) 

(775) 

(776) 

(777) 
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Further  calculate 


d2  =  D'  BS 2  tan*S 


HS2  =  B^j  tan 


S  2 


BR2  BS2  Sk 


DN2  =DN(D,d2,^s,Hs2,T,BR2,HR) 
and  use  Bs } ,  DNt ,  BS2  and  DN2  to  start  the  iterative  routine. 


(778) 

(779) 

(780) 

(781) 


Bs  i+2  Bs  i  +  dn. 


BS  i+  1  '  BS  i 


i+i  DN, 


(Nc  -  DN.) 


di+  2  ®  ”  BS  i+2  tan 


di+2  <  0  ? 


(782) 

(783) 

(784) 


If  yes,  execute  Equations  (785)  and  (786)  and  continue.  Otherwise,  ignore  these  two 
equations  and  continue. 


di+2  =  10  m 

(785) 

D-  10m 

Bs  I+  2  tan  ips 

(786) 

HS  i+2  =  BS  i+2  tan 

(787) 

BR  i+2  ~  BS  i+2 

(788) 

DNj+2  -  DN(D,  dj+2 ,  <ps ,  Hs j  +  2,  T,  BR  j+2 ,  HR ) 

(789) 

ei+2  =|Nc-DNj+2| 

(790) 

ej+2  <  0.00001  m  ? 

(791) 

If  yes,  stop  and  call  the  associated  dj+2  “d..”  Exit  from  the  iterative  routine.  Otherwise, 
reenter  the  iterative  routine. 
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End  of  the  iterative  routine. 


The  estimated  depth  of  the  seamount  peak  is 


d 


general  coin  |>. 


I)  -  BSj  tan  =  dj 


( 792 > 


where 


Jsj 


B 


S  i+2 


P93) 


Note  that  general  compensation  includes  the  total  absence  of  a  root  as  a  subcase.  To  avoid  singu¬ 
larities  in  Equations  (704)  and  (706),  one  should  not  attempt  to  realize  this  case  by  setting  HR 
to  zero  but  by  assigning  to  it  a  very  small  value  instead.  HR  =  0.000001  m  has  proven  to  be  a 
suitable  choice. 


Estimation  of  Peak  Depth  -  Output  Data 

For  the  isostatically  compensated  seamount,  store  and  print  out  (indicating  the  units) 


G/g 

m2/gr 

P  s 

gr/m3 

Pw 

gr/m3 

Pr 

gr/m3 

Pm 

gr/m3 

D 

m 

T 

m 

V 

PR 

,  deg 

W0 

km 
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isostat.  com  p . 


deg  and  deg-min-arc  sec 


deg  and  deg-min-arc  sec 


The  index  “j”  above  indicates  that  the  quantities  so  labeled  are  the  results  from  the  last  cycle  of 
the  iteration  performed  (values  associated  with  d.). 

For  the  generally  compensated  seamount,  store  and  print  out  (including  the  units) 


m2/gr 
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r 


Wo 

km 

Nc 

m 

dJ 

in 

”SJ 

in 

BSJ 

m 

brj 

in 

^general  comp. 

or  d  as  applicable 

uncom  p.  rr 

^SM 

deg  and 

deg-min-arc  sec 

^SM 

deg  and 

deg-min-arc  sec 

m 


The  index  “j”  above  indicates  that  the  quantities  so  labeled  are  the  results  from  the  last  cycle  of 
the  iteration  performed  (values  associated  with  d.). 


ANALYSIS  OF  DATA  DISPERSION 


Input  Data 


Obtain  the  following  from  the  SEASAT-A  data  track  and/or  data  track  header  (Tables  3 
and  4). 


oN  Calculated  average  of  all  oN  (Word  6,  Table  4)  on  data 

track  under  consideration.  This  is  the  average  geoid 
height  confidence  bound.  m 

Oj  Calculated  average  of  all  (Word  7,  Table  4)  on  data 

track  under  consideration.  This  is  the  average  vertical 
deflection  confidence  bound.  arc  sec 

vs  Satellite  subtrack  velocity  m/sec 

At  Time  interval  between  data  points  sec 
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Further,  have  available 


Re  Earth’s  “radius”  (Rt  =  6378000  m)  m 

G/g,  ’  P\v  *  Pr  *  Pm  ♦  as  specified  above, 

D  Ocean  depth  m 

AD  Uncertainty  in  ocean  depth  (user  input).  Normally,  use 

a  value  of  10%  of  D.  m 

T  Crustal  thickness  (user  input)  m 

AT  Uncertainty  in  crustal  thickness  (user  input)  m 

Seamount  slope  angle  deg 

Seamount  slope  angle  deg 

Nc  “Observed”  geoid  elevation  above  seamount  peak  m 

W0  Initial  estimate  of  seamount  width  as  base  m 

sk  Geometry  factor  for  generally  compensated  seamount 

cases 

Hr  “Height”  of  root  for  generally  compensated  cases  m 


Error  Estimates 


^lSM  *  At 

(794) 

ANSM  85  °N 

(795) 

A5  «  a5 

(796) 

AtA  *  At,  * 

At 

(797) 

AaA  555  AaB  » 

7r  A5  rad 

=  Ac* 

180  3600  arc  sec 

(798) 
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r 


deg 

A<fi.  A*p  =%!  (1.5E+05)  Aa - = 

a  t>  ratj  s 


(7Wl 


ana  ~  anb  =»  ansm  =»  on 


ANC  *=  -  0.3  Nc 


(800) 
(801  ) 


..  .  ^  vs  At  180 

axsm  51  a*sm  %  ~£ - ~dez 


(802) 


AW,.  *  2.8  vs  At 


(803) 


Data  Dispersion  for  lsostatically  Compensated  Seamounts 


Symbolize  the  computer  routine  that  performs  the  peak  depth  estimation  for  isostatically 
compensated  seamounts  by 


d. 


isostat.  comp. 


=  f(D,  T,  *pA 


*b>Wo 


Nc) 


(8U4) 


Let 


^isostat.  compjp  +  Ap 

symbolize  an  evaluation  of  the  above  function  (computer  routine)  f( . )  with  all  arguments 

nominal  except  the  parameter  “p”  indicated.  The  latter  shall  enter  the  routine  with  a  value  per¬ 
turbed  by  its  error  increment.  Calculate 

^iso st a t,  compj  D  +  A  D 
^isostat.  comp.^T  +  AT 

^isostat.  comp,^A  +  A^^  ^isostat.  comp.^g  +  A*?B 
^iso stat.  comp.^  Wq  +  AW^ 

^isostat.  compJ  +  AN^, 
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and 


(^d)  ^isostai.  comp  J  D  +  AD  disostat.  comp. 

(Ad)^j  ^isostat.  comp.^T  +  At  disostat.  comp. 

~  (disostat.  comp.^A  B  +  A^ia  „  -  disoslal  com 

(AU)*W0  *  co«p.)W0  +  AWq  -  disosUt  comp 

(Ad)‘iNc  *  (diso sta t.  eomp.)Nc  ♦  Anc  -  djsomt  comp 


Data  Dispersion  Analysis  for  Generally  Compensated  Seamounts 

Symbolize  the  computer  routine  that  performs  the  peak  depth  estimation  for  generally 
compensated  seamounts  by 


general  comp.  k(  D,  T,  ) 

and  calculate 

(Ad)^D  ^general  com  p  ?  D  +  AD  ^general  com  p. 

(^d)Ai  ^general  comp.^  T  +  At  ^general  comp. 

(^)A^  ^general  comp.^  ^  g  +  A^^  g  ^general  comp. 
^d)AW0  ^  ^general  co  m  p  )  W  Q  +  AW  Q  ~  ^general  com  p. 
(Ad>AN^  ^  ^general  comp.^  +  ANc  ”  ^general  comp. 


Additional  Details  Concerning  the  Data  Dispersion  Analysis 

Note  that  the  two  algorithms  indicated  by  Equations  (804)  and  (810)  depend  on  parameters 
other  than  the  ones  listed.  For  example,  the  densities  are  missing,  sk  and  HR  are  omitted  from 
Equation  (810).  All  these  quantities  were  disregarded  because  the  necessary  quantitative  data 
on  the  structure  of  the  seamount/root  system  and  on  seamount  composition  were  not  available 
with  the  degree  of  accuracy  that  would  make  inclusion  into  a  data  dispersion  study  meaningful. 
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AUTOMATIC  ANALYSIS  OF  SATELLITE  ALTIMETRY 


To  perform  an  automatic  seamount  survey  on  SEASAT-A  altimetry,  first  specify  the  desired 
data  track  per  device  set  number  and  revolution  number.  Be  aware  that  a  SEASAT-A  altimetry 
track  is  likely  to  consist  of  several  distinct  track  segments.  For  each  of  these  track  segments,  the 
present  seamount  detector  will  automatically  perforin  the  Matched  Filter,  the  Roughness  Detector, 
the  Seamount  Locator,  and  the  Calculation  of  Initial  Seamount  Parameters.  For  each  track 
segment,  indicate  the  cases  rejected  and  those  retained  for  further  processing.  For  the  cases  re¬ 
jected  as  well  as  for  those  retained,  print  out  the  longitude  XSM  and  the  latitude  y?SM  resulting 
from  Equations  (415)  and  (416). 

For  the  valid  detections,  execute  the  Estimation  of  Peak  Depth  for  Isostatically  Compensated 
Seamounts  plus  the  associated  Analysis  of  Data  Dispersion.  Then  execute  the  uncompensated 
subcase  (sk  =  2  and  HR  =  0.000001  m)  of  the  Estimation  of  Peak  Depth  for  Generally  Com¬ 
pensated  Seamounts  plus  the  associated  Analysis  of  Data  Dispersion. 

Make  provision  for  the  user  to  execute,  automatically  for  the  entire  data  track  or,  optionally, 
for  selected  individual  detections,  the  Estimation  of  Peak  Depth  for  Generally  Compensated 
Seamounts  and  the  associated  Analysis  of  Data  Dispersion  for  arbitrarily  specified  values  of 
sk  and  HR . 
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APPENDIX  A 


TRIAL  DATA  FOR  COMPUTER  PROGRAM  CHECKOUT 


Test  Trajectory:  SEASAT-A  Revolution  No.  1375 

This  is  a  North-East  to  South-West  pass  leading  over  Seamount  Gregg  in  the 
New  England  Seamount  Province  and  over  Bermuda. 

Test  Object:  Seamount  Gregg  in  the  New  England  Seamount  Province. 

Dimensions  of  Test  Object: 

ds  *  1400  m 

%  8.7  deg 

Ws  %  47  km 

XSM  *  61  deg  W  =  299  deg  E 

^sm  %  38  9  deSN 

Source:  Unclassified  “CONTOUR  SHEET  0707  (NAR-6)”  made  available  to 
the  author  in  July  1979  by  CODE  3513,  Naval  Oceanographic  Office,  NSTL 
Station,  Bay  St.  Louis,  Miss.  This  is  a  bathymetric  contour  map  containing 
Seamount  Gregg. 

Geophysical  Constants: 


G/g  = 

0.68024 

E-14 

m2/gr 

on 

II 

2.60 

E+06 

gr/m3 

Pw  = 

1.03 

E+06 

gr/m3 

Pr  = 

2.95 

E+06 

gr/m3 

Pm  = 

3.40 

E+06 

gr/m3 

Procedure  for  Checkout  Calculations: 

The  seamount  detector  computer  program  was  applied  to  SEASAT-A  Rev. 
No.  1375.  Seamount  Gregg  was  detected  in  due  course.  The  computer  ex¬ 
ercised  all  segments  of  the  above  described  algorithm.  In  particular,  estimations 
of  peak  depth  were  performed  for  the  isostatic  case,  the  uncompensated  case, 
and  a  generally  compensated  case  characterized  by  specified  values  for  sk  and 
Hr.  A  complete  data  dispersion  analysis  was  done  for  each  of  these  three 
cases. 
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To  verify  the  computer  results,  the  entire  algorithm  from  Program  Segment 
“Calculation  of  Initial  Seamount  Parameters”  on  was  coded  by  the  author 
for  use  on  a  programmable  electronic  calculator.  Starting  with  the  output  data 
of  Program  Segment  “Seamount  Locator,”  the  rather  vast  system  of  calculator 
routines  was  exercised.  All  computer  results  were  in  the  process  duplicated 
to  a  sufficient  number  of  digits. 


ERROR  ESTIMATES 


II 

0.041243731 

m 

®6  = 

0.83135406 

arc  sec 

ifiS  = 

6.7840 

km/sec 

At  = 

0.490024 

sec 

AD  = 

10%  of  D  =  0.1  D 

AT  = 

600 

m 

A*>s  *  A*a/b 


7r  deg 
1.8  x  3.6 


Of, 

- —  =0.40305  deg 

arc  sec 


9.3081  km 
0.4493  m 

TRIAL  DATA  FOR  ISOSTATIC  COMPENSATION 


AW0  *  2.8  vs  At 
ANC  w  -  0.3  Nc 


Nominal  Case: 

Data  prior  to  estimation  of  peak  depth: 
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(Ref:  Program  segment  input  data  plus  Equations  (720)  through  (729) 


D 

= 

5000  m 

T 

= 

5000  m 

Nc 

= 

1.4977448  m 

= 

9.8951328  deg 

(BS>0 

= 

20711.482  m 

(ds)o 

= 

1387.080712  m 

(Hg)0 

= 

3612.919288  m 

(Br)o 

= 

20711.48200  m 

(Hr>o 

= 

12605.07396  m 

(DN)0 

0.849684309  m 

Results  of  estimation  of  peak  depth: 

(Ref:  Equations  (754),  (750),  and  (753)) 

(B^  =  26493.53167  m 

(DN^  =  1.497744767  m 

(ds>j  -  <duomt.comp.>„ominal  =  378.4576320  m 

Perturbed  Cases: 

(For  each  case,  the  single  parameter  indicated  is  perturbed.) 

1.  AD  =  10%  of  D  =  500  m 

D  -*•  D  +  AD  =  5500  m 
(Bs)j  =  27105.91440  m 
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3. 


4. 


(DN^  =  1.497744736  m 

”  ^isost.  comp.^D  +  iD  ~  1 .6333  It i 

(Ads)AD  =  393.1757  m 
AT  =  600  m 
T  -*■  T  +  AT  =  5600  m 
( B s).  =  26157.90074  m 
(DN).  =  1.497744800  m 

<ds>j  =  (diSo,..comp.)x  +  AT  =  437.0052  m 

(Ads)AT  =  58.5476  m 

A^s  =  0.40305  deg 

spt  •*  sPs  +  A^>s  =  10.2981828  deg 

(Bs).  =  25771.9928  m 

(DN).  =  1.49774479  m 

(dj)j  ^diso »t.  com  p +  Ai^g  —  317.2795  m 
(AdSW  =  -61.1781  m 


AW0  =  9.3081  km 

(Bs)0  -*•  (Bs)0  +  ~ AWq  =  25365.5340  m 
(Bs).  =  26493.5321  m 
(DN).  =  1.49774483  m 


(ds)j  =  (d, 


is o st .  com  p.^fi 


SO 


W 


=  378.4576  m 


(Ads)Aw^  =0.0000  m 
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5. 


ANC  =  0.44932344  m 
Nc  ->  Nr  +  ANC  =  1.04842136  m 
(B^  =  22682.94818  m 
(DNlj  =  1.0482  m 


(tig  )j  ^isost.  comp.^N^,  +  AN^, 

(Ads)ANc  =  664.7198  m 


1043.1774  m 


TRJAL  DATA  FOR  UNCOMPENSATED  CASE 


Nominal  case: 

Data  prior  to  estimation  of  peak  depth: 

(Ref:  Program  segment  input  plus  Equations  (763)  through  (771)) 


D  =  5000  m 

T  =  5000  m 

Nc  =  1.4977448  m 

y>s  =  9.8951328  m 

(Bs)0  =  20711.482  m 

(d$)0  =  1387.080712  m 

(H$)0  =  3612.919288  m 
sk  =2 

(Br)q  =  41422.96400  m 

Hr  =  0.000001  m 

(DN)0  =  1.871038679  m 

45 


* 


Results  of  estimation  of  peak  depth: 


(Ref:  Equations  (793),  (789),  and  (792)) 

(B.)  =  18891.02732  m 

a  j 

(DN).  =  1.497744794  m 

(iL).  =  (d  )  .  .  =  1704.641562  m 

v  S  j  uncomp.nominal 

Perturbed  Cases: 

(For  each  case,  the  single  parameter  indicated  is  perturbed.) 

1.  AD  =  500  m 

D  -*■  D  +  AD  =  5500  m 

(Bo).  =  19219.5224  m 

(DN).  =  1.49774479  m 

(ds)j  =  (du„comp.)D  +  AD  “  2147.3388  m 

(Ads)iD  =  442.6972  m 

2.  AT  =  600  m 

T  ->  T  +  AT  =  5600  m 
(B„).  =  18891.0273  m 

“  J 

(DN)j  =  1.49774479  m 

«V>j  =  <duncon,p.>T+AT  =  1704.6416  m 

(Ads  )A  T  =  0.000  m  (Root  absent .  .  .  effect  of  crustal 
thickness  absent) 


r 


3.  A*>s  =  0.40305  m 

SPS  -^S+A^S  =  10.2981828  m 
(Bs  =  18539.6599  m 
(DN)j  =  1.49774479  m 

(ds)j  =  (duncomp.Vs  +  A^s  =  1631.3802  m 
t^s)^s  =  -73.2614  m 
4-  AW0  =  9.3081  km 

(Bs)0  ->  (B$)0  +  ~  AWfl  =  25365.4340  m 
(B^  =  18891.0273  m 
(DN)j  =  1.49774479  m 

(ds>j  =  (duncomp.»Bsn  +  -^w  =1704.6416  m 

o  u  j  0 

(Ads)Aw  q  =  0.0000  m 
5.  ANC  =  -0.4493  m 

Nc  ->  NC+ANC  =  1.04842136  m 
(Bs)j  =  16332.7632  m 
(DN)j  =  1.04842145  m 

<dS>j  =  (du„comp.)Nc+ANc=  3150.9062  m 
(Ads>ANc  =  446.2646  m 
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TRIAL  DATA  FOR  GENERALLY  COMPENSATED  CASE 


Nominal  Case: 

Data  prior  to  estimation  of  peak  depth: 

(Ref:  Program  segment  input  plus  Equations  (763)  through  (771)) 


D 

5000 

m 

T 

5000 

m 

Nc  = 

1.4977448 

m 

9.8951328 

deg 

20711.482 

m 

(ds)o  = 

1387.080712 

m 

(Hs)0  » 

3612.919288 

m 

sk  = 

2 

(Vo  = 

41422.96400 

m 

X 

73 

II 

3700  m 

(DN)0  = 

0.92643443 1 

m 

Results  of  estimation  of  peak  depth: 

(Ref:  Equations  (793),  (789),  and  (792)) 

(B-)  =  23952.54037  m 

2*  J 

(DN)j  =  1.497744778  m 

(d  )  =  (d  ,  )  .  ..  =  821.7091810 

'uS'j  '  general  com  p/nom  inal 
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Perturbed  Cases: 


(For  each  case,  the  single  parameter  indicated  is  perturbed.) 


AD 


500  m 


=  1251.5463  m 


D  -*■  D  +  AD  =  5500  m 


(B«).  =  24354.7574  m 
(DN).  =  1.49774473  m 


^gen.  comp.^D  +  AD 

(Ads)AD  =  429.8371  m 
AT  =  600  in 
T  ->T  +  AT  =  5600  m 
(B,,).  =  23849.9808  m 

^  J 


(DN).  =  1.49774479  m 


(dS)]  ^gen.  compA  +  AT  839.5997  m 

(Ads)AT  =  12.8905  m 

Ai£s  =  0.40305  deg 

t^s  -►  >fis  +  A<ps  =  10.2981828  deg 


(Bo ).  =  23339.7922  m 

a  3 

(DN).  =  1.49774480  m 


(ds)j  ^gen.  eomp.Vs  +  Ai^s  759.2056 

(Adg  =  “ 62.5036  m 


m 
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4. 


AW0  =  9.3081  km 


(Bs)0  -  (BS)0  +  -AW0  =  25365.5340  m 


(Bs)  =  23952.5403  m 


5. 


(DN).  =  1.49774477  m 


(d,)8  =  (d. 


S'j  '  gen,  comp.'BS0  q 

(Ads  )AWq  =  0.0000  m 
ANC  =  -0.4493  in 


)  *  w  =  821.7092  m 

’h„„  +  n 


Nc  -►  Nc  +  ANC  =  1.04842136  m 


(Bq  ).  =  21475.1599  m 

^  J 

(DN)j  =  1.048421461  m 


(dc).  =  (d  )M  .  Aw  =  1253.8644  m 

v  S  v  gen.  com  p/N^  +  AN^ 

(Ads)ANc  =  432.1552  in 


50 


DISTRIBUTION 


Defense  Mapping  Agency 
Naval  Observatory  Building  56 
ATTN:  O.  W.  Williams 
C.  Martin 
P.  M.  Schwimmer 
W. Senus 

R.  Berkowitz 
Washington,  DC  20360 

Defense  Mapping  Agency 
Aerospace  Center 
ATTN:  L.  B.  Decker 

M.  Schultz 
K.  Nelson 

St.  Louis,  MO  63118 

Defense  Mapping  Agency 
Hydrographic/Topographic  Center 
6500  Brooks  Lane 
ATTN:  H.  Heuerman 
Washington,  DC  20315 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  VA  22314 

Library  of  Congress 

ATTN:  Gift  and  Exchange  Division 

Washington,  DC  20540 

Commander 

Naval  Oceanographic  Office 
NSTL  Station 
ATTN:  T.  Davis 
J.  Hankins 

S.  Odenthal 

Bay  St.  Louis,  MS  39522  f 


DISTRIBUTION  (Continued) 


NASA  -  Wallops  Flight  Center 
ATTN:  J.  McGoogan 

Wallops  Island,  VA  23337 

NASA  -  Goddard  Space  Flight  Center 
ATTN:  J.  G.  Marsh 

R.  Kolenkiewicz 
Greenbelt,  MD  20771 

Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
ATTN:  W.  E.  Giberson 
P.  Rygh 

A.  Loomis 
J.  Lorell 
G.  Bom 

Pasadena,  CA  91103 

NOAA  -  National  Ocean  Survey/National  Geodetic  Survey 
6001  Executive  Blvd. 

ATTN:  B.  H.  Chovitz 

B.  Douglas 

C.  Goad 

Rockville,  MD  20852 

NOAA  -  Pacific  Marine  Environmental  Laboratory 
ATTN:  J.  Apel 
M.  Byrne 
Seattle,  WA  98105 

Aerospace  Corporation 
2350  East  El  Segundo  Boulevard 
ATTN:  Library 

El  Segundo,  CA  90245 

Applied  Physics  Laboratory 
Johns  Hopkins  University 
8621  Georgia  Avenue 
ATTN:  John  McArthur 
C.  Kilgus 

Silver  Spring,  MD  20910 


DISTRIBUTION  ^Continued) 


Oceanographer  of  the  Navy 
Naval  Observatory  Building  1 
Washington,  DC  20360 

Naval  Oceanography  Division 
Naval  Observatory  Building  1 
ATTN;  Code  NOP-952  (H.  Nicholson) 
Washington,  DC  20360 

Deputy  Chief  of  Naval  Operations 
ATTN:  NOP-21  ID 
Washington,  DC  20350 

Strategic  Systems  Project  Office 
ATTN:  SP20123  (P.  Fisher) 
Washington,  DC  20376 

Office  of  Naval  Research 
800  N.  Quincy  St. 

ATTN:  P.  C.  Badgley 
Arlington,  VA  22217 

Naval  Research  Laboratory 
ATTN:  V.  Noble 
P.  Vogt 

Washington,  DC  20375 
Commander 

Fleet  Numerical  Oceanography  Center 
Monterey,  CA  93940 

Department  of  the  Air  Force 
HQ  Space  &  Missile  Test  Center  (AFSC) 
ATTN:  Sven  C.  Ernberg 
Vandenberg  Air  Force  Base,  CA  93437 

NASA  Headquarters 

600  Independence  Ave.,  S.W. 

ATTN:  W.  F.  Townsend 
Washington,  DC  20546 


DISTRIBUTION  (Continued) 


University  of  Texas 
Dept,  of  Aerospace  Engineering 
ATTN.  B.  D.Tapley 
B.  Schultz 
Austin,  TX  78712 

Scripps  Institution  of  Oceanography 
ATTN:  R.  Stewart 

R.  L.  Bernstein 
LaJolla,  CA  92093 

•  Woods  Hole  Oceanographic  Institute 

ATTN:  J.  A.  Whitehead 
Woods  Hole,  MA  02543 

Institute  of  Marine  and  Atmospheric  Science 
City  University  of  New  York 
ATTN:  W.  J.  Pierson 
New  York,  NY  10031 

Department  of  Geodetic  Science 
Ohio  State  University 
ATTN:  R.  Rapp 

Columbus,  OH  43210 

Local: 

E41 

K05 

K10 

K104  (Ugincius) 

K104 (Zondek) 


K 1 2  (40) 

K12(Groeger)  (40) 

K1 2  (Smith) 

K1 1 
K13 

K14  (5) 

X210  (6) 


